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Abstract 

We investigate the scalar field perturbations of the 4+1-dimensional 
Schwarzschild black hole immersed in a Godel Universe, described by 
the Gimon-Hashimoto solution. This may model the influence of the 
possible rotation of the Universe upon the radiative processes near a 
black hole. In the regime when the scale parameter j of the Godel 
background is small, the oscillation frequency is linearly decreasing 
with j, while the damping time is increasing. The quasinormal modes 
are damping, implying stability of the Schwarzschild-Godel space-time 
against scalar field perturbations. The approximate analytical formula 
for large multipole numbers is found. 


1 Introduction 


Black hole’s behavior is often crucially dependent upon the cosmological 
background in which the black hole is immersed. The simple case of a 
black hole immersed in an asymptotically flat space-time is described by 
the Schwarzschild solution. A natural extension is to consider a cosmo¬ 
logical term which is described in terms of the Schwarzschild-de Sitter or 
Schwarzschild - anti-de Sitter solutions. Recent investigations of radiative 
processes around such black holes show that the radiative features are de¬ 
pendent upon the asymptotic conditions on infinity. Thus, for example, the 
quasinormal spectrum and late time tails are totally different for Schwar¬ 
zschild pQ, Schwarzschild-de Sitter and Schwarzschild - anti-de Sitter 
black holes |3]. 

The cosmological background we are interested in here, is the rotating 
Universe. The rotation seems to be a universal phenomenon: all compact ob¬ 
jects in the Universe rotate. Yet the standard Friedman-Robertson-Walker 
metric represents rather idealized model of isotropic homogeneous world 
filled with perfect fluid. It looks improbable that such a finely tuned universe 
can exist starting from the Big Bang to the present stage. In the beginning 
of the investigations of rotating cosmological models it was suggested that 
one should observe the anisotropy of the Microwave Background Radiation 
(MBR), yet, as shown later, that the rotating models with no anisotropy of 
MBR or broken causality can exist |^. In addition, apparent anisotropy in 
distribution of the observed angles between the polarization vectors and po¬ 
sition of the major axis of radio sources can be related to a possible rotation 
of the Universe [3]. For further advance in the possibility of observation of 
global rotation see the review [Hj. 

An exact solution for the rotating Universe was found by Godel [7]. 
His solution was originally proposed for a four dimensional space-time. It 
possesses, among others, the following properties: it is homogeneous, has 
rotational symmetry, and allows the definition of the direction of positive 
time consistently in the whole solution, and, what was in the focus of further 
research, it allows closed time-like curves, i.e. the time machine. 

Recently, the Godel Universes have been of considerable interest jHj, be¬ 
cause in five-dimensional minimal supergravity the maximally supersymmet¬ 
ric Godel-type universes are U-dual to pp-waves. Thereby, the Godel-type 
universes are important as a an opportunity of quantizing strings in this 
background a well as due to its relation to the corresponding limit of super- 
Yang Mills theories. On the gravitational side, the pp-waves dual to the 
Godel Universe corresponds to the Penrose limit of near-horizon geometries. 
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As far as we are aware, until recently, an exact solution for a stationary 
black hole immersed in a rotating Universe was not known. Nevertheless, 
such a solution has been obtained by Gimon and Hashimoto within the 
above mentioned five-dimensional minimal supergravity [HI- This solution 
represents the 4-|-1-dimensional Schwarzschild space-time when the scale pa¬ 
rameter of the Godel background j goes to zero, and to the five dimensional 
Godel Universe when the black hole mass vanishes, thereby giving us the 
model for the Schwarzschild black hole immersed in the rotating Universe. 
The different features of this solution have been investigated recently in a 
series of papers m- The generalizations of the Gimon-Hashimoto solution 
were obtained in m- 

Yet, here, we are interested in this solution from a rather different point 
of view, namely, we would like to find out what will happen with black hole 
(classical) radiation in the rotating Universe. The straightforward way to 
know it, is to investigate the quasinormal modes which govern the black 
hole response to external perturbations at late times. The quasinormal 
spectrum is sensitive to boundary conditions both at the event horizon and 
at spatial infinity, so the spectrum must be considerably affected by the 
rotating cosmological background. 

In the present paper we had to be limited by the case of “slow rota¬ 
tion”, i.e. the case when the influence of the cosmological background is 
weak. This happens when the above mentioned parameter j is small. We 
found that at least in the regime of small j, the black hole is stable against 
scalar field perturbations and as a result all found modes are damping. Due 
to cosmological rotation the real oscillation frequencies are decreasing and 
are roughly proportional to j, while the damping rates are decreasing non¬ 
linear ly with j. In addition we derive an approximate analytical formula for 
QN modes with large multipole number L. Fortunately the regime of small 
j seems to be the most reasonable phenomenologically. 

The paper is organized as follows: in Sec.2 we give some introductory 
material on the Schwarzschild-Godel metric. In Sec 3., the scalar field equa¬ 
tions is obtained in the limit of small scale parameter j of the cosmological 
background. For any j, the Klein-Gordon equation is not separable at least 
in the coordinates we have considered. In Sec.IV we find the quasinormal 
frequencies for the scalar field perturbations. In the end we discuss the 
obtained results and future perspectives. 
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2 Preliminaries of the Schwarzschild-Godel space- 
time. 


The bosonic fields of the minimal (4+1)- supergravity theory consist of the 
metric and the one-form gauge field, which are governed by the equations 
of motion 

R^u = 2{F^o.F^-^g^,F^'] ( 1 ) 


^ 2\/3 






( 2 ) 


Here, SaX-y/iu \/ 's^ax-y/iu- 

In the Euler coordinates {t,6','ip', cj)'), the solution of the equations of 
motion (1), (2), describing the Godel universe, has the form (Hj: 

= —{dt + + — (dd'^ + + dcj)'"^ + 2cos9'dil^'d4>'), (3) 

where (t£ = dcj)' + cosOdilj'. The parameter j defines the scale of the Godel 
background. At j = 0 we have the Minkowski space-time. The solution for 
the Schwarzschild black hole in the Godel universe is given by [2] 

3^2 


ds"^ = —f{r)dr — g{r)ra'ldt + h{r)r‘^ + k{r)dr^+ 

2 

^{d6'‘^ + dip''^ + dcj)'"^ + 2cos9' dY' d(j)'), 


(4) 


where 


f{r) = 1 - 


2M 


g(r) = 2jr, 


HY = + 2M), k{r) = ( 1 - + 


2M 16j2M2 


-1 


(5) 


The radius of the event horizon is also corrected by parameter j. 


rsH = x/2M{l — SpM). 


( 6 ) 


Note that the maximal value of the black hole mass M is l/8j^. For a larger 
mass the horizon area vanishes and one has a naked singularity. 

The above black hole metric keeps five of the nine isometries of the 
Godel universe, generated by dt, and by four generators of the SU (2) x U (1) 
subgroup of the SO {A) isometry group acting on [2j. 

In the limit j = 0 we have the (4+l)-dimensional Schwarzschild solution, 
while in the limit of m = 0 the pure Godel space-time is recovered. To treat 
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the scalar field perturbations around such a Schwarzschild-Godel black hole 
let us rewrite the metric in the bi-spherical coordinates {4>, ?/>, 9), which are 
connected with the Euler angles (cp', ip', 9') by 

4>'= Ip + 4>, ip' = ip-(p, ,9'= 29. (7) 

Then, in the regime of small j, i.e. discar ding terms of order the metric 

takes the form 


Ri —f{r)dt^ — 2g{r)r{{sin9)^d(p + {cos9)‘^d'ip)dt + k{r)dr^+ 
r‘^{d9‘^ + {cos9)‘^dpP‘ + {sin9)'^ dpP'). 

Note that in the above equation k{r) = (l — 

Up to the inverse metric has components 

-1 


11 / 2A^ 

d'' = - 1 - 


nn 2M oo 1 

9 =^-— - 


( 8 ) 


( 9 ) 


= r--\cos9)-\ g^^=r-\sin9)-\ g^-= g^^ =-2j [ I - 


-1 


( 10 ) 


3 Scalar field perturbations of the Schwarzschild - 
Godel space-time. 

The scalar held perturbations in a curved background are governed by the 
Klein-Gordon equation 

□ ( 11 ) 

Since the background metric has the Killing vectors dt, 9^, the wave 
function T can be represented in the form 

^ ^ ^i^t+ik^+im4>Y{^g)R{^r). ( 12 ) 

Unfortunately, variables in the Klein-Gordon equation are not separable, 
at least in the considered coordinates for the full Gimon-Hashimoto metric. 
The separability is connected with the existence of the Killing tensor na, 
and, it is possible that one can separate variables in some other privileged 
coordinate systems. Here we were limited to small values of j, for which the 


5 



variables in the Klein-Gordon equation can be separated even in ordinary 
bi-spherical or Euler coordinates. 

Using the expressions for metric coefficients (8-10), the scalar field equa¬ 
tion (11) takes the form: 


A 

dr 


1 ( 1 - 


2M\ ^dR{r) 


A 


.+{uj +4:ju){k+m))f ^{r)R{r) + ^R{r) = Q, 
or j 

(13) 

where the separation constant comes from the equation for angular variables, 


1 


cosOsinO 89 \ 


^ (sinOcosO^'] - 
86 




+ 


m 


{cos9)^ {sinOy 


Y = \Y. 


(14) 


Going over to the tortoise coordinate dr* = dr/f{r) and to the new 
wave function ik = i?(r)r^/^, the equation (13) can be reduced, after some 
algebra, to the wave-like form 


+(j? + 4:juj{k + m) -V (r*)^ vk = 0. 
The effective potential has the form 

V{r) = /(r) + ^fir) + 


f ^ I ^ I T + m){2l Y k Y m Y 2) 


(15) 


(16) 


Here I, k, and m run over the values 0,1, 2,... The tortoise coordinate r* is 
defined on the interval (—oo,-|-oo) in such a way, that the spatial infinity 
r = -|-oo corresponds to r* = oo, while the event horizon corresponds to r* = 
—oo. The above effective potential is positively defined and has the form 
of the potential barrier which approaches constant values at both spatial 
infinity and event horizon. In fact, the potential V (r) coincides with that 
for (4-|-l)-dimensional Schwarzschild black hole when taking the multipole 
number L to be 2Z -|- A: -|- m, yet, the spectrum is different, due to the term 
Ajuj{k Y m) in (5), which depend not only on the final value L, but also on 
terms I, k, m. 


4 Quasinormal modes of the Schwarzschild-Godel 
black hole. 

If choosing a positive sign for the real part of w (w = Reco — ilmoj), QNMs 
satisfy the following boundary conditions 

'k(r*) ~ Ci exp(±zt(;r*), r—> Too, (17) 
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corresponding to purely in-going waves at the event horizon and purely out¬ 
going waves at infinity. 

In order to find quasinormal frequencies of the black hole with an effective 
potential in the form of the potential barrier (16), we use the WKB approach. 
The WKB approach was used for calculations of quasinormal modes in the 
first order beyond eikonal approximation by Schutz and Will m, extended 
by Iyer and Will to the third order m, and recently extended to the sixth 
WKB order The WKB approach up to the 6 th WKB order has been 
used used recently in a series of papers m, m, [iHi, where QN frequencies 
of different black holes were considered, and, the comparison with accurate 
numerical values showed very good agreement. The accuracy of the WKB 
results is the better, the larger the multipole number L and the less the 
overtone number n. In fact for n larger then L the WKB formula cannot be 
applied. 

From here we shall use the units such that 2M = 1. 

The WKB formula has the form m- 

i ° - h 2 - T 3 - L 4 - Ls - Le = n -I- -, (18) 

v'— 2 Fq ^ 

where Vq is the height and Vq is the second derivative with respect to the 
tortoise coordinate of the potential at the maximum. L 2 and L 3 can be 
found in in, L 4 , L 5 and Lg are presented in m, the corrections depend on 
the value of the potential and higher derivatives of it at the maximum. 

The 6 th order WKB values of the quasinormal frequencies are shown 
on Figures 1-4 and Table I. From Fig. 1 one can see that the real part 
of uj is decreasing with growing of j, being roughly proportional to j, for 
a fixed black hole mass and fixed values of I, m, k. This can be easily 
explained in the following way: from the wave equation (15) one can learn 
that if one discards small values of order O(j^), then + Aj{m + k)u! = uJq, 
where u)q is the Schwarzschild value of u under some hxed M, /, m, k, n. 
Furthermore this can be represented as {u)+ 2j{k->rm))‘^ — 8 j^(/c-|-m)^ = ujq, 

i. e. the real part of iv is roughly increased by 2j{k + m), while the change in 
imaginary part comes from the term 8j'^{k + m)^. More accurately, Fig. 2 
shows that the imaginary part is decreasing when j is increasing. Therefore, 
the influence of rotating cosmological background, represented by parameter 

j, gives rise to decreasing of the oscillation frequency and of the damping 
rate. Thus, in the rotating Universe the QN modes damp more slowly, but, 
because of the considerable falling down of Reuj and slight falling down of the 
Imuj, the resulting quality factor Reoj/21 muo is decreasing and, thereby, the 
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black hole in a rotating Universe is a worse oscillator than in a non-rotating 
one. 

In Table 1. we put the low overtones of the QN spectrum for different 
values of k, m, I, and j. From that table we can learn that the higher any of 
the values k, I, or m and the lower the overtone n, the better the accuracy 
of the WKB formula, and, as a result, the less is the difference between 
the 6th and 3th order WKB data. An intrinsic fact for any black hole 
quasinormal spectrum, when the overtone number grows ReuJ falls down, 
while the damping rate grows. We cannot judge what will happen with 
asymptotically high overtones (n ^ oo), since we analyze here only an 
approximate solution. 

In the eikonal (high frequency) approximation, we can use the first order 
WKB formula for finding the lower overtones. Thus, for large I, and thereby 
for large L = 21 + k + m, in units 2(2M)“^ we obtain; 

L +1 ,, , 2n -|-1 , , 

UJ = ^— + 2j{k + m) -i . (19) 

Note that there are two limitations on this formula. First, when k or m is 
also large, the general relation 

9 X fL + l 2n + l\‘^ , , 

Lv -4:j{k + m)uj= - i j (20) 


holds. Moreover, we should be careful when interpreting this formula at 
asymptotically large L, since it uses an approximate metric and correc¬ 
tions in metric may produce different asymptotic values. Yet we expect it 
should be correct for moderately large values of L. Note that the above for¬ 
mulas are accurate enough even for not very large values of L, for instance, 
for L = 4 the relative error is about several percents. 

The massless scalar QNMs analysis can easily be extended to the massive 
case, in which one has the same wave-like equation but with the effective 
potential 




r -t- 


( 21 ) 

Yet we can use the above WKB formula only for small values of the field 
mass fi, since for large fx the effective potential has three turning points. 
The 6th order WKB frequencies are presented in Fig.3 and 4 as functions of 
fi. Thus, we see that the larger the held mass, the larger is the real part of 
Lv, and the smaller the imaginary part. In other words, the “massive” QN 



modes decay more slowly and have greater real frequency of oscillation. It is 
known that at asymptotically high overtones the mass of the field does not 
affect the QN modes m- Note that all the above features were observed 
for massive scalar field of the Schwartzshild m and Reissner-Nordstrom 
black holes m, m •Since these features were found also for massive Dirac 
field (see m and references therein), it is possible that they are generic for 
massive fields of arbitrary spin. 

TABLE I: WKB values for QNMs 

at fixed j = 1/8. 2M = 1. 


k 

m 

1 

n 

3th WKB order 

6 th WKB order 

0 

0 

0 

0 

0.49123 - 0.41100ii 

0.54633 - 0.36087i 

1 

0 

0 

0 

0.78080 - 0.35384i 

0.79151 - 0.35575i 

1 

0 

0 

1 

0.59794 - 1.15113i 

0.61846 - 1.14902i 

0 

0 

1 

0 

1.50707 - 0.35787i 

1.51050 - 0.35770i 

0 

0 

1 

1 

1.38524 - 1.10754ii 

1.39249 - 1.10537i 

0 

0 

1 

2 

1.19442 - 1.90690ii 

1.18639 - 1.94829i 

1 

0 

1 

0 

1.77155 - 0.35330i 

1.77293 - 0.35317i 

1 

0 

1 

1 

1.67618 - 1.07942i 

1.67927 - 1.07838i 

1 

0 

1 

2 

1.51426 - 1.84424i 

1.50511 - 1.86038i 

1 

0 

1 

3 

1.30766 - 2.63920i 

1.27639 - 2.73256i 

1 

1 

1 

0 

2.05408 - 0.34832i 

2.05475 - 0.34826i 

1 

1 

1 

1 

1.97396 - 1.0582K 

1.97552 - 1.05774i 

1 

1 

1 

2 

1.83110 - 1.79890i 

1.82449 - 1.80621i 

1 

1 

1 

3 

1.64511 - 2.5705K 

1.61726 - 2.61762i 

1 

1 

1 

4 

1.42547 - 3.36434ii 

1.37379 -3.51096i 


5 Discussions 

We have investigated the decay of (generally speaking, massive) scalar field 
around a Schwartzshild black hole immersed in a rotating cosmological back¬ 
ground. In the limit of the small cosmological parameter j, the QNMs, which 
govern the decay of the scalar field at late times, have been found. It was 
found that the cosmological rotation gives rise the decreasing of the real 
frequencies of oscillations (proportional to the cosmological parameter) and 
of damping rates. The quality factor of the black hole as an oscillator is 
smaller in the presence of cosmological rotation. The massive scalar field 
damps more slowly and have greater oscillation frequency. All found modes 
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are damping what supports the stability of the Schwartzshild-Godel space- 
time against scalar held perturbations. Yet, within the approximate solution 
we analyzed, one cannot judge about stability eventually. Note also that the 
stability of the metric as such is determined by the gravitational perturba¬ 
tions, although the scalar held perturbations may coincide with tensor type 
gravitational perturbations m which are decisive in gravitational stability 
m The present analysis can also be extended to the case of scalar held 
interacting electromagnetically with the charge of the black hole uni , i.e. to 
the case of the decay of charged scalar held around a Reissner-Nordstrem- 
Godel black hole. 
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Figure 1: Real part of u; as a function of the Godel background scale pa¬ 
rameter j (diamond k = 2, m = 1, I = 1), (box A: = 1, m = 1, / = 1), (star 
k = 0, m = 1, I = 1). 


Im (oj) 



Figure 2: Imaginary part of u; as a function of the Godel background scale 
parameter j (star k = 2, m = 1, I = 1), (box k = 1, m = 1, I = 1), (diamond 
k = 0, m = 1, I = 1). 





Figure 3; Imaginary part of w as a function of the mass jj, for ^ = 1, A: = 1 
m = 1 (box), I = 1, k = 0, m = 1 (star), I = 1, k = 0, m = 0 (diamond) 

j = 1/8. 
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Figure 4: Real part of w as a function of the mass n for I = 1, k = 1, m = 
(box), / = 1, A: = 0, m = 1 (star), 1 = 1, k = 0, m = 0 (diamond); j = 1/8. 







